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Many features of real granular fluids under rapid flow are exhibited as well by a system 
of smooth hard spheres with inelastic collisions. For such a system, it is tempting to 
apply standard methods of kinetic theory and hydrodynamics to calculate properties 
of interest. The domain of validity for such methods is a priori uncertain due to the 
inelasticity, but recent systematic studies continue to support the utility of kinetic theory 
and hydrodynamics as both qualitative and quantitative descriptions for many physical 
states. The basis for kinetic theory and hydrodynamic descriptions is discussed briefly 
for the special case of a low density gas. 



1. Introduction 

Granular media in rapid, dilute flow exhibit a surprising similarity to ordinary fluids 
and the utility of a hydrodynamic description for such conditions has been recog- 
nized for many years [Q. This phenomenology has come under scrutiny in recent 
years with questions about the domain of its validity and the associated consti- 
tutive equations appearing in the hydrodynamic equations H,^. Answers to such 
questions can be found in a more fundamental microscopic description where the 
tools of nonequilibrium statistical mechanics are available for critical analysis. An 
intermediate mesoscopic description between statistical mechanics and hydrody- 
namics is that of kinetic theory, whose applicability to granular matter also poses 
questions. The objective in this short communication is to provide an example of 
a precise context in which the relevance of kinetic and hydrodynamic descriptions 
can be assessed. A more general review of the current status for this problem with 
extensive references was prepared for this Workshop and should be consulted 
for completeness. 

The system considered is a one component gas of N smooth hard spheres at low 
density. The inelastic collisions are characterized by a normal restitution coefficient 
a < 1, where a = 1 is the elastic limit. The state of the system at each time 
is specified by a point in the 6N dimensional phase space and the dynamics is 
given by uniform motion and binary inelastic velocity changes at contact. These 
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are the necessary ingredients to construct a statistical mechanics for this idealized 
model of granular media. It is straightforward to write the Liouville equation for 
the probability density and from it obtain the BBGKY hierarchy for the reduced 
distribution functions • In the next section the condition of low density is 

exploited to obtain a complete, closed kinetic theory description for the gas from 
the BBGKY hierarchy. The low density expansion leading to this result is formally 
the same as that for a gas with elastic collisions. Consequently, it would appear 
that the kinetic description for the granular gas has the same level of validity as 
that for normal gases. 

The derivation of hydrodynamics from the kinetic theory is discussed in 
the third section. Hydrodynamics is defined generally as the composition of exact 
balance equations for the density, energy density (or granular temperature), and 
the momentum density (or flow field) plus constitutive equations for the associated 
fluxes and cooling rate. Constitutive equations exist whenever there is a "normal" 
solution to the kinetic theory. Such a solution can be constructed explicitly for 
weakly inhomogeneous states, leading to a Navier-Stokes hydrodynamics for the 
granular gas |]9|, p^JTT| . The need to go beyond the Navier-Stokes approximation is 
rare for normal gases, but is not uncommon for many relevant granular gas states. It 
remains a challenge for granular gas theory to understand the form of constitutive 
equations for these more general conditions, but this uncertainty should not be 
interpreted as a signature of the failure for hydrodynamics to apply. 



2. Kinetic Theory 

Concerns about a kinetic theory description for granular gases are often based on 
a too restrictive concept of the prerequisites (e.g., the existence of an equilibrium 
state) . It is useful therefore to provide a formal " derivation" of the kinetic descrip- 
tion at low density to demonstate the close similarity between normal and granular 
gases, without unwarranted conceptual prejudices. In this section a small param- 
eter expansion of the reduced distribution functions is shown to provide a formal 
solution to the entire BBGKY hierarchy, in parallel to the corresponding analysis 
for normal gases | [T^ . The analysis is similar to an expansion proposed by Grad for 
the hard sphere gas with elastic collisions jl^. In this approach, there is no ref- 
erence to concepts such as " approach to equilibrium" , Maxwellian distribution, or 
"molecular chaos", and the distinction between inelastic or elastic collisions plays 
no explicit role. Thus, superficially at least, it appears the basis for the low density 
kinetic theory is the same in both cases. 

The s-particle reduced distribution functions, /^*^(a;i, • • ■ ,Xs,t), obey the BBGKY 
hierarchy where Xi — (q,, v^). A dimensionless form of this hierarchy is obtained by 
scaling the space and time with the mean free path £ = l/(na'^) and the mean free 
time to = £/vf). Here, n is the density, a is the hard sphere diameter, and wq is some 
characteristic velocity. Similarly, the reduced distribution functions are scaled with 
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(u/vq)^ . The resulting diniensionless BBGKY hierarchy has the form 

^ s ' _ \ 

i=l i<j J 

= J dxs+iTii,s + l)f^'+'\x,r--,Xs+i,t). (2.1) 

The binary scattering operator T{i^j) for a pair of particles {i^j} is defined by 
T(i, j) = J da e(gy • a)ig,j ■ a) [a-^S{q,j - Aa)6,r.i _ ^^q^^. ^ ^^^^ ^ ^2.2) 

where tr is a unit vector along q,;j = q-i — and da denotes a two dimensional 
solid angle integration over the sphere for particles at contact. Also, = Vj — Vj 
is the relative velocity, and is the scattering operator defined for any function 
X{-Vi,Wj) by 6~ X(vi,Vj) = X(v-,v^). The "restituting" velocities are 

/ l + a, , 1 + a, ... 

V,; EE V, - ^—(gy •cr)cr, Vj- EEVj + -^(gij •cr)cr. (2.3) 

The a dependence of T{i,j) contains all aspects of the inelasticity, and plays no 
explicit role in the following expansion. 

In this dimensionless form the BBGKY hierarchy depends on the single dimen- 
sionless parameter \ = a/£ = na^, the ratio of the "force range" to the mean free 
path. This parameter is small at low density, suggesting an expansion for a solution 
to the entire hierarchy as a power series in A. The dependence on A occurs explic- 
itly as shown on the left side of ( ^.iD and implicitly through the finite separation 
of the colliding particles in T(i,j). The structural features of the expansion in A 
are simplest if it is performed at fixed T{i,j). The s-particle reduced distribution 
functions are taken to have the representation 

/(^) , • • • , a;,, t) = /^^^ (a:i , • • • , x„ t) + A^/f ^ (a^i , • • • , x„ i) + . .. (2.4) 

It is then readily shown that the hierarchy is solved exactly to order in the form 

s 

f^''\xi,---,xs,t)^]lf^'\x,,t), (2.5) 

i=l 

s s s s 

fi'\x^r ■ ■ ,xs,t) =J2l[ f!,'\x„t)f['\x,,t)+Y, n fo'\^k,t)G{x,,x„t),{2.6) 

(s) 

where the expression for holds for s > 2. Thus, the reduced distribution func- 
tions for any number of particles are determined as a sum of products of the single 
particle functions fQ^\xi,t) and f[^\xi,t), and pair function G(a;i, a;2, i). These 
are determined from the set of three fundamental kinetic equations 

(I + VI • Vi) f^'Hxut) = J(xi,i I /(^^ (2.7) 
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I- + VI . Vi - /i + V2 • V2 - /2 ) G{XI , X2 , t) = T(l, 2)f^^^ (Xi , i)/^'^ {X2 , t) (2., 



dt 



d_ 

dt 



Vi-Vi-/i^/(')(xi,t) = J dx2Til,2)G{xi,X2,t), (2.9) 



Here J{xi,t \ /o^"*) is the Boltzmann-Bogoliubov collision operator and /i, defined 
over functions of xi, is its linearized form 

Jixi,t I /(^)) = y dx2T{l,2)f^'\xi,t)f^'\x2,t) (2.10) 
Jift(xi)= f dx2T{l,2)(^f^^\xi,t)h{x2) + hixi)f^^\x2,t)y (2.11) 



These low density results (2.7)-(2.9) provide the kinetic theory description for 



the gas. They are remarkably rich. As expected, the leading order distribution func- 
tion /q^^ is the solution to the Boltzmann equation. The two particle correlations 
are generated from the uncorrelated product of Boltzmann solutions through inelas- 
tic binary collisions T{l,2)fQ^\xi,t)fQ^\x2,t) on the right side of (|2^). Finally, 
corrections to the Boltzmann solution due to correlations are given by a coupling 
of the distribution function to the correlations in ( |2.9| ) (the so-called "ring" rec- 
oUision effects). The solutions and implications of these kinetic equations can be 
quite different for elastic and inelastic collisions. But these differences come from 
the equations themselves and should not be interpreted as signatures of their fail- 
ure to apply. For example, at a = 1 a possible solution for an isolated system 
is f^^\xi,t) /Af(wi), G{xi,X2,t) = = f[^\xi,t), where Jm is the Maxwell- 
Boltzmann distribution. Equation (|2.9| ) supports fm because energy is conserved, 
and G{xi,X2,t) = because T(l, 2)/ji/(vi)/m(v2) = for the same reason. Since 
energy conservation no longer holds with a < 1 it is not surprising that the isolated 
system does not approach equilibrium, the Maxwellian is not a stationary solution, 
and that finite correlations exist. Indeed, the extent to which such predicted dif- 
ferences agree with observations from molecular dynamics provide support for the 
kinetic theory, not limitations on it as is sometimes implied. 

Clearly, the above derivation has not restricted this kinetic description to iso- 
lated systems or near-equilibrium states. In fact, the most interesting cases of prac- 
tical interest are response to boundary conditions and/or external fields. The sim- 
ilarities between normal and granular fluids is closest for such " nonequilibrium" 
conditions. Too often, properties of granular gases are contrasted only to those of 
the equilibrium state for normal gases. 

It is important to note that practical access to the solutions to the above kinetic 
equations is possible for a wide range of conditions by direct simulation Monte Carlo 
(DSMC) . Much attention has been given to the special " homogeneous cooling 
state" (HCS) which is a solution to the Boltzmann equation with the scaling form 

/^^)(v,<) - VoHt)ncj,{v/vo{t)) , vlit) = 2T{t)/m. (2.12) 
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The form of this is HCS distribution is known to good approximation by analytic 
methods, which have been confirmed by DSMC [||. The correlations for this state 
also have been studied in some detail [T5|] by both analytic and simulation methods. 
Finally, studies of more complex states (e.g., shear flow) also have been given JT6| . 
In summary, the kinetic theory description appears to describe well a wealth of new 
phenomena peculiar to inelastic collisions. 



3. Hydrodynamics 

Consider now a spatially inhomogeneous state, created either by initial preparation 
or by boundary conditions. The local balance equations for the density n(r, i), 
granular temperature T{r,t) (or energy density), and flow velocity U(r,t) follow 



directly by taking moments of the Boltzmann equation (2.7) with respect to 1, v, 
and 

A^^ + ?^V■U = 0, (3.1) 
DtT + ^ {P^J^,U, + V • q) = -TC, (3.2) 

DtU, + {mn)-^djP,, ^0. (3.3) 

Here Z3( = 9( + U • V is the material derivative, Pij{v, t) is the pressure tensor and 
q(r, t) is the heat flux. The form of these balance equations is the same as for fluids 
with elastic collisions except for the source term on the on the right side of ( ^^ ) 
due to the dissipative collisions, where C oc (l — a^) is identifled as the cooling rate. 
The fluxes Pij , q and the cooling rate C are given as explicit low degree moments 
of the distribution function f^-^\xi,t) 

P., =y"dvmy,y,/(r,v,i), q = J d^,^mV^Vf{r,^r,t), (3.4) 

C-(l-a')Y^ / ^virfv2d<7e(a-g)(a.g)3 /(ri,vi,i,)/(r2,V2,i). (3.5) 

The utility of these balance equations is limited without further specification of 
Pij, q, and C which, in general, have a complex space and time dependence. How- 
ever, for a fiuid with elastic collisions this dependence "simplifies" on sufficiently 
large space and time scales where it is given entirely through a functional depen- 
dence on the fields n, T, and U. The resulting functional dependencies of Pij and q 
on these fields are called constitutive equations and their discovery can be a difficult 
many-body problem. The above balance equations, together with the constitutive 
equations, become a closed set of equations for n, T, and U called hydrodynamic 
equations. This is the most general and abstract notion of hydrodynamics, which en- 
compasses both the Navier-Stokes form for small spatial gradients and more general 
forms for nonlinear rheological transport. The primary feature of a hydrodynamic 
description is the reduction of the description from many microscopic degrees of 
freedom to a set of equations for only five local fields. 
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The critical problem for a hydrodynamic description is therefore to determine 



the existence and form of the constitutive equations. It is clear from (3.2) and (3.4) 
that they can be obtained if the Boltzmann equation admits a "normal" solution, 
whose space and time dependence occurs entirely through its functional dependence 
on the fields 

f{r,v,t)^F{v\n,T,U). (3.6) 

The fluxes and cooling rate then inherit this space and time dependence and become 
constitutive equations. The space and time dependence of the fields follows from 
solution to the resulting hydrodynamic equations to complete the self-consistent 
description of F. An example is given by the HCS distribution in ( ^.12| ) where 
there is no space dependence and all of the time dependence occurs through T{t). 



The latter is determined from the hydrodynamic equation (3.2), which reduces to 



dtT = —T(. Use of ( 2.12 ) in (3.5) gives the constitutive equation C — CoT^^^ where 
Co is a constant. 

For gases with elastic collisions the prototypical hydrodynamics is that of the 
Navier-Stokes equations. The corresponding constitutive equations, Newton's vis- 
cosity law and Fourier's heat law, follow from a normal solution to the Boltzmann 
equation obtained from an expansion in the spatial gradients. The reference state 
is the local Maxwellian whose parameters are the exact density, temperature, and 
flow velocity for the nonequilibrium state being described. Deviations from this ref- 
erence state are proportional to the spatial gradients of the temperature and flow 
field. The systematic expansion for the normal solution to the Boltzmann equa- 
tion is generated by the Chapman-Enskog method. There have been three main 
objections/reservations/concerns regarding application of this method for granular 
gases: 1) the absence of an equilibrium state as the basis for the local Maxwellian 
reference state, 2) the inherent time dependence of any reference state due to colli- 
sional cooling, and 3) the inclusion of the energy (temperature) as a hydrodynamic 
field when it is not associated with a conserved density and does not have a time 
scale solely characterized by the degree of spatial inhomogeneity. 

The first two concerns are primarily technical rather than conceptual issues that 
can be answered by direct application of the Chapman-Enskog method to see if it 
indeed generates a normal solution to the granular Boltzmann equation. Consider 
a state for which the spatial variations of n, T, and U are small on the scale of the 
mean free path, ^Vlnn << 1, where £ = l/na'^ is the mean free path. Then it is 
expected that the functional dependence of the normal solution on the hydrody- 
namic fields can be made local in space through a Taylor series expansion about a 
point r for which the distribution function is being evaluated. Let e denote a formal 
small "uniformity" parameter measuring the small spatial gradients in the fields. 
It is worth emphasizing that there is no restriction on the cooling rate C; only the 
spatial gradients are being ordered by the uniformity parameter. The distribution 
function, collision operator, and time derivatives are given by the representations 

F = FW+eF« + ..., Js = + 9* = + e9« + ... (3.7) 
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The coefficients in the time derivative expansion are identified from the balance 
equations with a similar expansion for Py , q, and C generated through their defi- 



nitions (3.4) and ( |3.5| ) as functionals of F. As for elastic collisions the leading term 
F'"^ is further constrained to have the same moments with respect to l,w^, and v 
as the full distribution F. This assures that the hydrodynamic fields occuring in 
F^"^ are exactly those for the nonequilibrium state, and hence F^^"^ is normal. To 
zeroth order in e the Boltzmann kinetic equation becomes 

ic^"Vy (vfW) - jW[F("',fW]. (3.8) 

where V = v — U. 

The first concern above regarding the reference distribution now can be ad- 
dressed. The distribution F*^"^ is not free to be chosen, but rather is determined by 
the kinetic equation itself. For the granular gas it is not the local Maxwellian, but 



rather the local (normal) HCS solution (i.e., (|2.12| ) with the density, temperature, 
and flow field replaced by their nonequlibrium values). There is no a priori require- 
ment of an equilibrium state for the Chapman-Enskog method to apply, and in fact 
early applications based on a local equilibrium state are inconsistent and lead to 
incorrect transport coefficients even for weak dissipation. The concept of " approach 
to equilibrium" is no longer relevant for granular gases. 

Mathematically, the changes in this method for granular gases arise from the 
fact that the time derivative of the temperature does not vanish to lowest order in 
e, as it does for a gas with elastic collisions. The reference state F^^^^ incorporates 
this zeroth order time dependence of the temperature even for strong dissipation. 
This is the origin of the second concern noted above. However, since F^''^ is normal, 
it necessarily has the exact time dependence of all hydrodynamic fields. The only 
difference for granular gases is the introduction of a new time scale 1/^'-*''' in the 
reference state. There is nothing a priori inconsistent with a description of slow 
spatial decay towards a time dependent reference state. In fact, this is an important 
feature of the Chapman-Enskog scheme for both granular and normal fluids alike - 
use of a time independent reference state would limit the derivation to only linear 
hydrodynamics. 

Implementation of the Chapman-Enskog method to the flrst order in e is now 
straightforward and has been carried out in detail and without approximation re- 
cently for the Boltzmann equation ^ (and for its dense fluid generalization, the 
Enskog equation ||l^); the case of a two component mixture is considered in [ pT| . 
The constitutive equations for the one component fluid found to this order are 

P^J ^ p6^J ' Ti {d,U, + d,Uj - ■ - V • U, (3.9) 

q^-KVT-^Vn, (3.10) 

The form of the pressure tensor is the same as that for fluids with elastic collisions, 
where rj (a) is the shear viscosity depending on the restitution coefficient. The 
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heat flux is similar to Fourier's law, where k (a) is the thermal conductivity. There 
is a new transport coeSicient fi (a) coupling the heat flux to a density gradient, 
which vanishes at a = 1. The transport coeflicients are given in terms of solutions 
to inhomogeneous integral equations involving the linearized Boltzmann collision 
operator. Solubility conditions for the existence of solutions have been proven and 
approximate solutions in terms of polynomial expansions have been obtained for 
practical purposes, just as for the case of elastic collisions ||l^. Consequently, the 
transport coefflcients are known as explicit functions of a. The Chapman-Enskog 
method places no explicit restriction on a so the results are not limited to weak 
dissipation. 

It remains to discuss the third concern, use of the temperature as a hydro- 
dynamic field which is no longer associated with a conserved density. Analysis of 
the Navier-Stokes approximation described above shows there are two classes of 
inherent time scales to this hydrodynammic description. The first scales with the 
wavelength of the phenomena considered and can be made large (long time scales) 
by considering sufficiently smooth disturbances. This is the case for gases with 
elastic collisions. For granular gases there is a second time scale in the temperature 
equation, the inverse cooling rate, which is bounded for a given value of the resti- 
tution coefficient. The concern is that this new time scale may be shorter than that 
required for validity of the normal solution to the kinetic equation. 

Qualitatively, a wide class of solutions for spatially inhomogeneous states evolves 
in two stages. During a short transient period of the order of the mean free time (the 
kinetic stage), the velocity distribution approaches closely the local distribution F^°^ 



and becomes normal as in (3^). Subsequently, on a longer time scale the space and 
time dependence of the distribution occurs only through the fields that are governed 
by hydrodynamic equations. This notion of a two stage relaxation is similar to that 
for gases with elastic collisions and is confirmed for granular gases by DSMC. The 
question, therefore, is whether the inverse cooling rate is smaller than the mean 
free time. Certainly, this is the case for weak dissipation since C^^-* is proportional 
to (l — Q^) . For strong dissipation it is necessary to study in detail the spectrum of 
non hydrodynamic (kinetic) modes to find the slowest mode and compare it with 
the inverse cooling rate. This has been done for the diverse examples of Brownian 
motion, kinetic models, and uniform shear flow (see Q). In all cases it is found 
that there is a separation of time scales between the slowest kinetic mode and the 
inverse cooling rate, even at strong dissipation. Based on these examples, it appears 
that inclusion of the temperature in the hydrodynamic description is justified for a 
wide range of degrees of dissipation. 

The validity of Navier-Stokes hydrodynamics and the dependence of the trans- 
port coefficients on the restitution coefficient has been verified in a number of simu- 
lations, both DSMC and MD, giving excellent agreement with the predictions from 
the Chapman-Enskog method. The tests at low density based on the Boltzmann 
equation have been reviewed recently by Brey and Cubero | p^ . 
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The Navier-Stokes equations equations associated with (3.9) and (3.10) are also 
known as Newtonian hydrodynamics. Although the Chapman-Enskog method can 
be carried out to second order in e (Burnett order), it is likely that failure of the 
Navier-Stokes approximation signals more complex non Newtonian behavior for 
which other methods to construct the normal solution are required that are not 
based on a small gradient expansion. 



4. Discussion 

Two questions have been addressed here: 1) can kinetic theory provide a valid 
mesoscopic description of rapid flow (fluidized) granular media?, and 2) can a hy- 
drodynamic description be formulated and justified for a macroscopic description? 
The origin of the kinetic description for a low density gas has been shown to follow 
from an expansion that makes no reference as to whether the collisions are elastic 
or inelastic. The latter affects only the solution to these equations. It appears that 
many of the questions, concerns, and objections raised regarding the validity of 
a kinetic theory description are removed in this way. Furthermore, kinetic theory 
provides a powerful tool for analysis and predictions of rapid flow gas dynamics 
when combined with the method of DSMC for numerical solution. Its full potential 
for both conceptual and practical questions has not yet been exploited. 

An important use of the kinetic theory is to determine if a hydrodynamic de- 



scription applies and to define its domain of validity. Equations (2.6) - (2.8) leave 
no room for speculation. For given initial and boundary conditions the solution 
either approaches a normal form on some space and time scale or it does not. The 
issue of a hydrodynamic description is now a precise mathematical question. This is 
particularly important if the hydrodynamic description is more complex than that 
of Navier-Stokes. Non Newtonian hydrodynamics is rare or unphysical for simple 
atomic systems with elastic collisions. In contrast, they it is more common for 
granular gases in steady states where the gradients are strongly correlated to the 
coefficient of restitution. Kinetic theory is perhaps the only systematic means to 
determine the constitutive equations in these cases. 

It is clear from recent studies that granular media exhibit a wide range of in- 
teresting phenomena for which a Navier-Stokes hydrodynamics is an accurate and 
practical tool (see Q for specific examples and references). Although the context 
here has been limited to low density gases a similar kinetic theory basis has been 
developed for Navier-Stokes hydrodynamics in moderately dense fluids [|o|. More 
generally, non Newtonian behavior puts granular fluids in a class of complex ma- 
terials with mysterious, as yet unexplained, properties Kinetic theory and 
hydrodynamics (in the broader sense) can be expected to provide much of this 
explanation. 
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